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Abstract. We present synchronization and contractivity estimates for the kinetic Ku- 
ramoto model obtained from the Kuramoto phase model in the mean-field limit. For 
identical Kuramoto oscillators, we present an admissible class of initial data leading to 
time-asymptotic complete synchronization, that is, all measure valued solutions converge 
to the traveling Dirac measure concentrated on the initial averaged phase. If two initial 
Radon measures have the same natural frequency density function and strength of cou- 
pling, we show that the Wasserstein p-distance between corresponding measure valued 
solutions is exponentially decreasing in time. This contraction principle is more general 
than previous L^-contraction properties of the Kuramoto phase model. 



1. Introduction 

The objective of this paper is to present a contraction property of the kinetic Kuramoto 
equation in the Wasserstein metric. The synchronization phenomena exhibited by various 
biological systems are ubiquitous in nature, e.g., the flashing of fireflies, chorusing of crickets, 
synchronous firing of cardiac pacemakers, and metabolic synchrony in yeast cell suspensions 
(see for instance [H S] ) • Winfree and Kuramoto |19l [30] pioneered the mathematical treat- 
ment of these synchronized phenomena. They introduced phase models for large weakly 
coupled oscillator systems, and showed that the synchronized behavior of complex biological 
systems can emerge from the competing mechanisms of intrinsic randomness and sinusoidal 
couplings. The kinetic Kuramoto equation has been widely used in the literature [1] to 
analyze the phase transition from a completely disordered state to a partially ordered state 
as the coupling strength increases from zero. Suppose that g = g{i^) is an integrable steady 
probability density function for natural frequencies with a compact support (see (|2.ip for 
details). Let / = f{9,0,,t) be the probability density function of Kuramoto oscillators in 
9 £ T := ]R/(27rZ) with a natural frequency at time t as in [20]. The kinetic Kuramoto 
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equation(KKE) is given as follows: 

dtf + de{uj[f]f) = 0, {e,n)eTxR, t>0, 

u[f]{9,n,t) = Q- K [ sm{e -e^)p{e^,t)de,, p{e^,t):= [ fdn*, ^^'^^ 

subject to the initial data: 

fi6,n,o) = foie,n), f hde = g{n). (1.2) 

JT 

Note that KKE (jl.ip can be regarded as a scalar conservation law with a nonlocal flux, 
and it has been derived from a previously proposed Kuramoto model [U [20]. However, to 
the best of authors' knowledge, few studies have investigated the qualitative properties of 
the KKE, such as an asymptotic behavior and stability of some equilibria. 

The main results of this paper can be summarized as follows. First, we present sufficient 
conditions for the emergence of completely synchronized states. More precisely, when many 
coupled limit-cycle oscillators have the same natural frequency (identical oscillators), and 
the support of the initial Radon measure is confined in a half circle, we show that the 
measure-valued solution approaches a multiple of the Dirac delta measure concentrated on 
the initial average phase value and the common natural frequency value time-asymptotically; 
in other words, complete phase synchronization occurs asymptotically. For this purpose, we 
lift the finite-dimensional result for the Kuramoto model (KM) to the infinite-dimensional 
KKE. Second, we present a contraction property of the KKE in the Wasserstein p-distance 
for measure valued solutions with the same natural frequency distribution by using a strat- 
egy similar to the one described in [TJ [21]. We define a cumulative distribution function 
of a density function / for the KKE, say -F, and we derive a new integro-differential equa- 
tion using its pseudo-inverse function. Then, we use simple techniques for the optimal mass 
transport in one-dimension, i.e., the equivalence relation between the Wasserstein p-distance 
and the L^-distance of the corresponding pseudo-inverse of F in order to obtain the expo- 
nential decay estimate of the Wasserstein p-distance between two measure-valued solutions. 

The rest of this paper is organized as follows. In Section 2, we briefly review the Kuramoto 
model and its mean-field version (the KKE), and we provide several a priori estimates. In 
Section 3, we revisit an existence theory of measure valued solutions to the KKE, and we 
present several a priori estimates, in particular, we provide a finite-time stability estimate 
for measure-valued solutions in a bounded Lipschitz distance. In Section 4, we show the 
emergence of completely synchronized states by lifting the corresponding results for the 
KM to the KKE using the argument of the particle-in-cell method. This strategy has been 
employed in the Cucker-Smale fiocking model in (6]. Section 5 is devoted to the contraction 
property of the KKE using the method of optimal mass transport as described in [Tj [21] . 
Finally, the conclusions are stated in Section 6. 
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2. Preliminaries 

In this section, we briefly review the particle Kuramoto model and its kinetic mean- 
field model. Consider an ensemble of sinusoidally coupled nonlinear oscillators that can be 
visualized as active rotors on the circle S^. Throughout the paper, we will identify a rotor 
with an oscillator. Let Xi = e^^^^ be the position of the i-th rotor. Then, the dynamics 
of Xi is completely determined by that of phase 6i. In the absence of coupling, the phase 
equation for 9i is simply given by the decoupled ODE system: 

JO 

-} = ni, i.e., Oiit) = 9i{o) + n,t, 

at 

where fij is the natural phase- velocity (frequency) and is assumed to be a random variable 
extracted from the density function g = g{ft): 

g{—0,) = g{0,), spt((7) is bounded 

(2.1) 



/ ng{n)dn = o, / g{n)dn = i. 



In the seminal work [19] of Kuramoto, he derived a coupled phase model heuristically from 
the complex Ginzburg-Landau system. The KM is given by 

dff K ^ 

subject to initial data: 

em = 9m. (2.3) 
Note that the first term on the R.H.S. of (j2.2p represents the intrinsic randomness, whereas 
the second term denotes the nonlinear attractive coupling. Hence, synchronized states for 
system (j2.2p will emerge, when the nonlinear coupling dominates the intrinsic randomness. 

The system (|2.2p has been extensively studied over the last three decades, and it remains a 
popular subject in nonlinear dynamics and statistical physics (see review articles and a book 
[DElIiailTlESlEH]). In[l8l[l2], Kuramoto first observed that in the mean-field limit (iV 
oo), the system (12. 2p with a unimodal distribution function g{i^) (which is assumed to 
be one-humped and symmetric with respect to mean frequency := J2i=i ^0 ^ 
continuous dynamical phase transition at a critical value of the coupling strength Kcr '■ 

2 

Kcr = — r^, in the mean-field limit. 
ngiO) 

Moreover, he introduced an asymptotic order parameter r°° € [0, 1] to measure the degree 
of the phase synchronization in mean-field limit: 

r'^(K) := lim lim — Ve^^'^*) , 

i=l 

and he observed that this quantity r°° changes from zero to a non-zero value, when the cou- 
pling strength K exceeds a critical value Kcr- Note that for an initial phase configuration 
that is uniformly distributed on the interval [0, 27r), the quantity r°° is exactly zero, whereas 
for a completely synchronized configuration 9i = 9c for i = 1, - ■ ■ ,N, r°° becomes the unity. 
Therefore we can regard r°° as the order parameter'^ measuring the degree of synchro- 
nization. Before we conclude this subsection, we recall a complete phase synchronization 
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result from [13]. For this purpose, we introduce the diameters of the phase configuration 
6 = (9i, ■ ■ ■ , 9i\i) and the natural frequency set as follows. 

Deit) := max \9i{t) - 9j{t)\, t>0, Dq := max |0j - %|. 

i<i,j<N i,j 

Proposition 2.1. [13j Suppose that the natural frequencies, the coupling strength and initial 
configuration satisfy 

Qi = iy^j, K>0, Do := De{0) < tt, 

and let 9 = 6{t) be the smooth solution to the system ()2.2p -( f2r3]l with initial phase configu- 
ration 9q. Then we have 

e-^*Do < De{t) < e^^"*I)o, t > 0, (2.4) 

where a is the positive constant only depending on the diameter of the initial phase config- 
uration given by 

sin Dq 
a := ——. 
Do 

We also recall the estimate of existence of a trapping region for non-identical oscillators 
from [9] as follows. 

Lemma 2.1. [9] Let 9 = 9{t) be the global smooth solution to (I2.2p - (l2.3p satisfying 



Then we have 



< Do < TT, Dn>0, K>Ke:-- . ^ . 

sm Vq 



[i] sup Dg{t) < Do <Tr. 

[ii) 3 to > such that s\xpDg{t) < D^ 

t>to 



where D°° is defined by 



D°° := arcsin 



Dn 



K 

Remark 2.1. // we set the average phase and natural frequency of the particles as 

N N 
1=1 i=l 

then from the particle KM ()2.2p . one can easily obtain 

9P{t) = 9P{0) + nPt, for all t > 0. 

Without loss of generality, we may assume that 0? = using the phase-shift framework. 
Then we notice that identical and non-identical oscillators satisfying the assumptions in 
Proposition \2.1\ and Lemma \2.1\ satisfy 

{9i{t) ^c(O) as t — 7> oo, for identical oscillators, 
9i(t) e (61^(0) - D°°, 9P{0) + D°°) for all t > to, for non-identical oscillators. 



Note that the conditions and decay estimates (|2.4p are independent of the particle-number N . 
For the related synchronization estimates for the KM, we refer to [9l [lOl [TTl [HI [161 ISIJ [23] . 
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We rewrite the system (|2.2p as a dynamical system on the extended phase space T x M 
for {9i,Qi) : For i = I,-- - ,N, 



Throughout this paper, we will use the interval [0,2-k) to denote T = ]R/27rZ, i.e., 9 € 
[0, 27r) implies that 9 satisfies 9 + 27rZ = 9. 



In this section, we briefly review the existence of measure valued solutions to (jl.ip . For 
the KM, the rigorous mean-fleld limit was first done by Lancellotti [20] using Neunzert's 
general theory for the Vlasov equation |24l I27j . Optimal transport arguments allow to 
generalize these results in several ways for granular and fiocking models [11 [5]. H. Chiba 
recently obtained the same mean- field limit based on functional tools [8]. For a later use 
and reader's convenience, we present several estimates for measure valued solutions to the 



3.1. Measure-theoretic framework. In this subsection, we discuss a measure-theoretic 
formulation of the KKE. 

Let A^([0, 2tt) x M) be the set of nonnegative Radon measures on [0, 2tt) x M, which can be 
regarded as nonnegative bounded linear functionals on C([0, 27r) x M). For a Radon measure 
€ -A^([0, 2tt) X M), we use the standard duality relation: 



The definition of a measure- valued solution to equation (jl.ip is given as follows. 

Definition 3.1. For T G [0,oo), let fi G L°°([0, T); A4([0, 27r) x R)) be a measure valued 
solution to (jl.ip with an initial Radon measure G A^([0,27r) xM) if and only if satisfies 
the following conditions: 

(1) is weakly continuous: 




t > 0. 



(2.5) 



3. Existence theory of measure valued solutions 



KKE. 




{^ti h) is continuous as a function of t, V /i G Co([0, 27r) x M). 
(2) fi satisfies the integral equation: V /i G CQ([0,27r) x M x [0,T)), 




(3.1) 



where uj[fi]{9,Q, s) is defined by 

uj[fi]{9,n,s) 



Q — K{fis * sin)6' . 



(3.2) 



Here * denotes the standard convolution, i.e., 
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Remark 3.1. (i) Let f = f{6,^},t) be a classical solution to (jl.ip . Then the measure 
fi := fdQdO is a measure valued solution to (ll.ip . 

(a) Note that the empirical measure 

N 



K = Ti^^0^(t) ^^n,{t), where (6'j(t), is a solution of (g^I) 



,N _ 

N 



i=l 

is a measure valued solution to (jl.ip in the sense of Definition 3.1. Thus, the solutions to 
the KM (j2.2p can he treated as measure valued solutions via an empirical measure. Here, 
5z^ is the Dirac measure concentrated at z = z^,. 

(Hi) Since the density function g{^) has a compact support and the dynamics (jl.ip only 
governs the 6-variable, we can see that the projected Q-support of iit also has a compact 
support as well. Under the assumption of the compact support of g, we can expand the class 
of test functions to C([0, 27r) x M). 

(iv) By choosing h = 0, in (13. ip (see above comment in (Hi)), we have 

Lemma 3.1. Suppose that the density function g = g{i^) has a compact support and the 
initial measure satisfies 

iflQ,^) = 0, 

and let ^ G L°°([0, T); A4([0, 27r) x M)) be a measure valued solution to (jl.ip . Then for 
t > 0, we have 

(^t,i) = (mo,i> = 1, {fit,e) = {fio,e), t>o. 

Proof. In (j3.ip . we set h = 1. Then, the R.H.S. of (j3.ip will be zero hence, we have 
conservation of total mass. For the time evolution of the first moment of 9, it follows from 
Remark 3.1 (iv) that 

{nt,n) = o, t>0. 

We now set h{e) = 6 in and use ([32]) to get 

{l-H,9) = {f^o,9)+ {i^Ls,uj[fis])ds 

t 



= {fio,9)-K {ij.s,{fis*sin)e)ds = {fio,9), 
Jo 

where we used the anti-symmetry of sm(9 — 9^) to determine that (/i^, {^s * sin)0) = 0. □ 

3.2. A priori local stability estimate. In this part, we recall the stability estimate for 
measure valued solutions to (II. ip in the bounded Lipschitz distance. This stability estimate 
is crucial to the global existence of a measure valued function for the KKE. First, we review 
the definition of the bounded Lipschitz distance presented in [151 EH EZ] . We define the 
admissible set <S of test functions as 



S 



|/i : [0,27r) X M^M : ||/i||loo < 1, Lip(/i) < l| , 



where 

up{h) ■= sup \hiOi,^i)-h{e2,n2)\ 
(6*1, 01)^(6*2,^2) K^ii^i) - {(^2,^2)1 

Definition 3.2. [Ml 127] Let fi,iy £ A4([0, 27r) x M) be two Radon measures. Then the 
bounded Lipschitz distance d{fj,, v) between fi and v is given by 

d{ii, u) := sup (/i, h) — {u, h) 

Remark 3.2. 1. The space of Radon measures A^([0,27r) x M) equipped with the metric 
d{-,-) is a complete metric space. 

2. The bounded Lipschitz distance d for compactly supported probability measures is equiv- 
alent to the Wasserstein-1 distance (Kantorovich- Rubinstein distance) Wi (see |27| j; 

Wi{ii,u):= inf r f r [ \{e-9,,n-n,)\^{d{9,n),d{e,,n,)), 

where is the set of all product measures on ([0, 27r) x R) x ([0,27r) x M) such that 

their marginals are n and z^, respectively. Both equivalent distances endow the weak-* con- 
vergence of measures with metric structure in bounded sets. 

3. For any h £ C([0,27r) x M) with ^ o.n'd Lip{h) < b, we have 

{n,h) — {iy,h) < max{a,b}d{fi,i'). 

3.3. Existence of a measure valued function for KKE. We briefly present the exis- 
tence of a measure valued solution to the KKE using a Lancellotti's argument |20j . We first 
note a simple approximation argument of measures by smooth positive densities. For given 
fiQ G L°°(0,r;>J([0,27r) x E)), we set for < e < 1, 

fil := ((1 - e)Mo + exR) *Ve= f " 0) Vei^ " ^) ((1 " e)Mo + exn){de, dQ), 

J[0,2tt]xR 

where % = irj^^rf^)., r]^ is a periodic compactly supported mollifier with period 2tt and 77^ 
is a standard compactly supported mollifier satisfying H^/e ||Li{[o,27r]) = 1) ll^ellLi{R) = 
and xu is the uniform probability measure on rectangle TZ enclosing spt(/Uo) such that 
TZ C [0, 27r] X [— C, C] with C > 0. It is straightforward to check that 

d{lj.f),l^o) '^Wi{fLl,fio) ^ 0, as e^O 

and spt(/2§) C [0, 27r] x [-0^,0^]. 

We remark that if spt(//o) C (0, 27r) x [— C, C], then spt(/iQ) c (0, 27r) x [—Ce,Ce] using 
a standard compactly supported mollifier r/g for sufficiently small e. Since /Iq is absolutely 
continuous with respect to Lebesgue measure dOdil. and with connected support, then for 
all e the initial approximated measure /Ig can be approximated by a Dirac comb of uni- 
form masses. This means that there exist point distributions {(^^q^, f^io^)}i=i,...,Af) whose 
dependence on is elapsed for clarity, such that 

°° i=i ' ' 



8 



CARRILLO, CHOI, HA, KANG, AND KIM 



Then we solve the KM with A^-particles: 

^r^-^E--(^^'^-^n ^ = 0' ^>0' ^ = l>--->^> (3-4) 



dt ' TV ^ ^ * J dt 



with initial data: 



With solutions {9^'^ (t) , Q^'^ (t)) of ()3.4p . the approximate solution /i^''^ for the measure 
valued solution can be constructed as a sum of Dirac measures, i.e., 

The results of [2^ ^^ \Tl\ imply the continuous dependence with respect to initial data 
for measure valued solutions. Therefore, there exists a constant C depending on T and the 
initial support of the measure /^O'/^o such that 

d{fil,fi^t)<Cd{fil,fil). (3.6) 

for all t G [0,T]. As a consequence of (j3.6p . we have the convergence of approximate 
solutions and particle approximations to the measure valued solution of (II. ip : 

dil^ulA'^) < C d{fLl,ii^'^) and d{jil,nt) < C d{flf^, fio). 
This obviously yields 

diilf^, fit) < C{d{ilQ'^, ill) + d{fio, flD) . 
Hence we first let e — )■ and then letting — )• oo to have 

d{flf^,fJ't) 0. 
Therefore we have the following theorem. 

Theorem 3.1. Suppose that the initial measure fiQ G A^([0,27r) x M) and let fi^ be the ap- 
proximate solution constructed by p.Sp . Then there exists a unique measure valued solution 
Ht G A^([0,27r) X R) to (jl.ip such that nt is the weak-* limit of the approximate solutions 
fif as N ^ oo, i.e., 

lim d{fit,fJ'f) = 0. 

From now on, let us assume that the initial measure is a smooth absolutely continuous 
measure with respect to Lebesgue with connected support. These assumptions can be 
eliminated by standard mollifier approximation as above. Therefore, we will proceed by 
working on smooth solutions and obtaining estimates depending only on quantities that 
pass to the limit in the weak-* sense, and thus, stable estimates under this approximation. 
To avoid too much repetition, this procedure will not be specified in the proofs below and 
the statements of the results will be written directly for measure valued solutions. 



4. Asymptotic complete synchronization estimate 

In this section, we present an asymptotic synchronization estimate for the KKE (ll.ip 
by Hfting corresponding results for the KM (j2.2p using the argument of the particle-in-ceh 
method [25] discussed in the previous section. 

Let fi G L°°([0,r);A^([0,27r) x M)) be a measure valued solution to ([TTT]) . and let R{t) 
and P{t) be the orthogonal 9 and ri-projections of spt(/_i() respectively, i.e., 

R{t) := Fesptint), P{t) := Pcspt(/it), 
Then it is easy to see that 

P(t) = P(0), t > 0. 

We also set 

Deifit) ■■= diam(ii(t)), Dn{fit) := diam(P(t)), M{t) := {fit, 1), 

We observe from Lemma |3. II that 

M{t) = {fit, 1) = (^0, 1) = M{0) = 1, 

and since r2c(0) = 0, we obtain 

ndt) = and 9c{t) = 9c{0), t > 0. 

Lemma 4.1. Suppose that the oscillators are identical and let fiQ he a given initial proba- 
bility measure in A1([0, 27r) x M) satisfying 

{f^o,9) = TT, De{fto)<TT, K>0. 

Then the measure valued solution fi to (jl.ip - ()1.2p with an initial datum fiQ satisfies 

Deifit) < De{fio)e-'^^\ t > 0, 

where 

_ _ sin Deifip) 

Proof. We use the approximation argument in Theorem 13.11 giving a initial particle approx- 
imation defined as in (|3.3p . Then, it follows from Proposition 12.11 that the approximate 
measure valued solution fit^ G 7V4([0, 27r) x R) satisfies 

Deiftf) < De{fi^)e-''^-\ t > 0, 

where 

sinDeifio) 

Since Theorem 13.11 implies that d{fit,fit^) — )• as N ^ oo. Hence, DQ{fi^) — )• DQ{fit) as 
— 7- oo and we obtain the desired result. □ 

Remark 4.1. Throughout the paper, without loss of generality, we assume that {fio,9) = vr 
in order to avoid any possible confusion arising from the periodicity of 9. In fact, if the 
oscillators satisfy the assumption in Lemma \4-1\ (or Lemma \5.1\) . the orthogonal 9 -projection 
of spt(fit), R{t) is confined to the interval (0, 27r) for all t > (see Remark \2.1\) . This 
property will also he significantly used in Section 5 (see Lemma [5^) . 



10 CARRILLO, CHOI, HA, KANG, AND KIM 

We now show that the measure valued solution to the system (jl.ip for identical oscillators 
will converge to a multiple of the Dirac measure concentrated on (^c(O), ^^c(0)) in the phase 
space {9, Jl). We set 

Theorem 4.1. Suppose that the oscillators are identical, and let /zq € A1([0, 27r) x M) he a 
given initial probability measure satisfying 

6c{0) = vr, Dg{fio) < TT and K > 0. 
Then the measure valued solution fit to (jl.ip with initial datum fj,Q satisfies 

lim d{jit,fJ'oo) = 0, 
where d = d{-, •) is the bounded Lipschitz distance defined in Section 3. 
Proof. Let h = h{9) G C([0, 27r)) be any test function satisfying 

||/i||L- < 1, Lip(/i) < 1. 
Then h can also be regarded as a test function in C([0, 27r) x M). Then we have 

h{9)fit{d9,dn)- [ h{e)noo{de,dn) 

[0,27r]xK y[0,27r]xI8 

\{e)flt{dd) - h{7r) < r \e - 7:\fit{de) < Dg{ii^)e-'^^K 
Jo 

where we used Lemma |4.H and fitidO) is a ^-marginal of the measure fit, i-e., 

MdO) := [ fit{de,dn). 

This implies that 

Kat 



d{fit, /^oo) < Dg{fio)e ^0, as t ^ oo. 



5. Stability estimate of the KKE 



□ 



In this section, we present the strict contractivity of measure valued solutions to the KKE 
by using the method of optimal mass transport [TJ \2T\ [29] . The strict contractivity result 
generalizes the L^-contraction result for the KM in [9]. 

5.1. Alternative formulation of the KKE. In this part, we derive an alternative form 
of the KKE, which is more convenient for deriving estimates in terms of the Wasserstein- 
distance. First, we study the existence of an invariant set for the KKE. 

Lemma 5.1. Suppose that the initial probability measure fiQ and the coupling strength K 
satisfy 

< Deifi^) < vr, < I?n(^o) < oo, K> .^"^^"^ . ■ 

sm De[fi\i) 

Then, there exist to > and G (0, ^) such that the measure valued solution /i to (II. ip 
with initial datum fiQ satisfies 

Deifit)<D°°, t>to. 
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Proof. We apply the argument similar to that in the proof of Lemma 14.11 Let > be 
given. Then we have the following approximation fi^ for /ig: 

1 ^ 

1=1 

We now solve the Cauchy problem for KM: 

d9 K 

— ^ = 0. 

dt 

subject to initial data (0^(0), Oj(0)) = {9iQ,Q,iQ). Theorem 13.11 implies that 

/if^) — )> as ^ oo , 

and thus, Dg{fi^) — > Dg{fit) and DQ{fi^) — > Dfi{fit) as — > oo. Hence we can take 
large enough such that Dq{^q) and Dq{ij,q) satisfies the conditions of Lemma [2TT1 Thus, 
we find that there exist > and D°°'^ such that 

Deifif) < D°°'^, t > t^, for A^ large enough, 

where 

- KsinDe{^^^)-Dn{^.^y - ^ K ^ ^ 2^ 

We now let A^ — )> oo to obtain the desired result. □ 

In the remainder of this section, from Remark 14.11 we assume that 

C (0, 27r) and t>0. (5.1) 

Under this assumption the solution is given by a smooth particle density function f{6, 0, t) 
m for all t > 0. For a given Q, we consider a one-particle density function / as a function 
of 6. Then we define the pseudo cumulative distribution function of /: 



F{e,n,t):=[ f{e^,n,t)de, {e,n,t) e[o,2Tr) x 

Jo 



and a pseudo-inverse cj) of F(-, Q,, t) function of 9: 

(l){Tl,n,t) := inf{0 : F{6,n,t) > ??}, rj G [0,g{n)]. 

As long as there is no confusion, we use the notation F~^[r], 0.,t) = (j) as the pseudo inverse 
of F as ^-function. Then it is easy to see that 

=7?. (5.2) 

Lemma 5.2. Let ^ he a measure-valued solution to p.ip - (|1.2p . and let cj) be the pseudo- 
inverse function of the cumulative distribution function F. Then we have 

(i) max{6l I 61 G i?(t)} = max (j){g{9),Q.,t). 

nespt{g) 

(ii) mm{9 \ 9 G R{t)} = min 0(0, J], t). 

Q.&spt{g) 

{Hi) max (j){g{Vt),Q.,t) - min 0(0, 17, t) < t > to- 

nespt{g) Q£Spt{g) 
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Proof. Since the estimate for (ii) is similar to that of (i) and the estimate for (iii) follows 
from the estimates (i) and (ii) , we only provide the proof for the estimate (i) . For notational 
simplicity, we set 

Om := max{e \ 6 € R{t)}. 
Then, by definition of /it, we have 

9m = max{9 \ 9 G sptg(/(6', f], t)) and Q € spt(5)}, 
where sptg(/(^, 0, i)) is the ^-projection of spt{f{9,Q,,t)). This yields 

9m = max{(/)(g(r2), $7, t) such that Q G spt((/)}, 
by definition of the pseudo-inverse function. This completes the proof. □ 

Next, we derive an integro-differential equation for the pseudo inverse <j). It follows from 
(fO) that the smooth solution f{9,n,t) to (fTTIl- ifO]) satisfies 



f{o,n,t) = 0, QeR, t>0. 
We differentiate the relation (|5.2p in t and use dgF = f to get 



dtFi9,n,t) +fi9,n,t) dt(i){ij,n,t) = Q. 



This yields 



1 



e=(i}{ri,n,t) 



e=(p{ri,n,t) 



f{9,n,t) 

n + K / sm{9^ - (I){r],n,t))f{9^,n^,t)d9^dn^ using dSl]) 
Jr Jo 

r rai^*) 

fl + K / sin((/)(r7*, $7*, i) — 0(7?, $7, t))(i7?*(i$7*. 



where we used 9^ = (l){r]^ , Q,^ , t) and relation ()5.2p to see f{9^:,Q^,t)d9^: = dr]^. Hence, the 
pseudo-inverse (p satisfies the following integro-differential equation: 

dt(j) = n + K / sin((/)^, — (j))dr]^dn^. (5.3) 



where we used abbreviated notations: 

(f>^ := <f>{r]^,n^,t), (f> := (f>{r],n,t). 

The following results is a simple consequence of the change of variables and Lemma 13.11 

Lemma 5.3. Let fit be a measure valued solution to - (|1.2|) with an associated pseudo- 
inverse function (j). Then, we have 

/ / (t)dridn= / / 9fit{d9,dn), -r (l)dT]dn = 0. 

Jr Jo Jr Jo Jr. Jo 
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5.2. Strict contractivity in the Wasserstein distance. In this part, we present the 
proof of the strict contraction property of the KKE. 

For the one-dimensional case, it is weh known [71 [29] that the Wasserstein p-distance 
Wp{fii, 1x2) between two measures fii and fi2 is equivalent to the L^-distance between the 
corresponding pseudo-inverse functions (pi and 02 respectively. Thus, we set 

Wp{fli,fl2){^,t) := \\ct)i{-,n,t) - 02(-,^^,OllLP{O,g{Q))i 1 < p < 00. 

Since Wp{fii, fj,2) depends on we introduce a modified metric on the phase-space {9,Q): 

Wp{i^ii,fi2){t) := \\Wp{ni,H2){-,t)\\LP(R), 1 < p < 00. 

Below, we assume that the density function g{il.) has compact support. Then, it is easy to 
see that Wp{ni, 112) is a metric that satisfies 



lim Wp{fii,fi2){t) = Woo(mi,^2)W, 

p— >-oo 



t > 0. 



(5.4) 



Recall that the sgn function is defined by 

sgn(x) = 



1, X > 0, 
0, x = 0, 
-1, x<0. 



Lemma 5.4. Let ^ be a measurable function defined on [0,g{Q)] x M satisfying 



TT 



|$(ry, n)\ < ^ and 



/ / <^>{7],n)dr]dn = 0. 
Jr Jo 



Then for 1 < p < 00, we have 











\^r],n)\P-^sgnmr],n)) - |$(?7,, f],)r is5n(^>(r/„ 



X sm( )ari^drjd\^a\l < 

2 vr 



Proof. For notational simplicity, we set 

^ := ^{T],n), := ^{r]^,n^), and 
A(r?,r?,,0,0,) := r|$|P-isgn($) - |f-^sgn($. 



■am 



\^{r])\Pdr]dn. 



sm 



and we decompose the domain [0, g{^)] x M as the disjoint union of three subsets: 
V:={{rj,n)m7],n)>0}, Z:={{7],n)\^{r,,n)=0}, AA:={(r/,0) I «I>(7?,J7) <0}. 

Then it follows from the condition Jq^^^ ^drjdO. = that 



I \^dr]dVL = I \^dr]dn 



(5.5) 
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We use [0, g{n)] xIR = Pu2:u7\Ato obtain 




JO 



H / +/ +/ +/ ]/S.{ri,ri^,VL,VL^)drjdr]^dQ.d^^ 

VxZ JMxV JvxV JAfxAf JzxZ 



distinct signs same signs 

We now consider the following sub-integrals separately. 

I{A,B):= A{T],'i]^,n,n^)dr]dr]^dndn^, A,Be{V,Z,M}. 

JaxB 

We claim the following: 



Case 


A 


B 


I{A,B) < 


I 


V 


Z 


J^\<^^\Pdri^dn^ 


II 


M 


z 




III 


z 


V 


-£^j^\^\pd7]dn 


IV 


z 


M 




V 


V 


M 


_1 

TV 

+iM\'^r'd 


^(^) Xv- \'^\^dT]dn + c{M) |$*|pd?7*df]* 

i]dn \^^\dr]^dn^ + \^*\P~^dr]^dn^ Jj^ \^\di]dn 




VI 


M 


V 


_1 

TT 

+ \<^\p~^dl 


C{M) \^vdr]dVt + C{V) y \^^\vdr]^dVL^ 
'jdn Jj^ \^^\dT]^dn^ + Jj^ \^^\f-'^dr]^dn^ \^\dr]dn 




VII 


V 


V 


_1 

TT 

-ivi^r'd 


^'P) Iv \'^\^dridn + C{r) \<^^\Pdr]^dn^ 
7]dQ fj, dr/^,(iri* — /p <I>* p'^^dr/*dr2* Jp \^\drjdQ, 




VIII 


M 


M 


_1 
TT 

-ij^m'-'d^ 


c{M) Jj^ \^\Pdr]dn + c{Af) Xv- \^*\Pd'n,dn, 

•]dn Jj^ \^*\dr]^dn^ - Jj^ \ ^^\P~'^di]^dn^ Jj^ \^\dr]dn 




IX 


z 


z 






where jC{A) denotes the Lebesgue measure of the set A: 

C{A) := / Idrjdn. 
J A 



We also note that 



C{V) + C{Z) + C{M) = / / ldr]dVt = / g{VL)dn = 1. 
Case I: In this case, we use the definition of A(r/, r/*, $7, 0*) and the inequality 



sinx > —X, for x G 

TT 



°-2, 



to determine that 



A(77, 77,, 17,^7,) = -|$,|P-isin^ < 

Z TT 
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This yields 

I{V,Z) <-- [ l^^lPdrjdrj^dndn^ = -^^^ [ 
^ JvxZ 71" Jv 

Case II - Case IV: The estimates are basicahy the same as in Case I. Hence, we omit 
their estimates. 

Case V: In this case, we have 

A{rj,r],,n,Q,) = + sin (^^^J^^^ 

< -^(i^r + i^*r + i$r>*i + i^*r>i). 

This yields the desired result. 

Case VI: Once we interchange V < — t- J\f, the same estimate holds. 

Case VII: In this case, we need to consider two subcases: 

Either <I) > > or 0. 
By considering each case, we have 

A(77,7y,,J],f],) = (|$|P-i-|(I>,|P-i)sin(^;l^) 

< - (\^r' - \<^,r') {\<p,\ - m 

vr 

This yields the desired result. 

Case VIII:: The estimate is exactly the same as in Case VII. Hence we omit its estimate. 
Case IX:: The estimate is trivial. 



We now add all cases and use (|5.5p to find 




A{ri, 77*, Q, ft^,)dridr]^dQdQ^, 
< (C{V) + C{Z) + C{M)) / \^\Pdr]dn 



vr 







\^\PdT]dn. 

□ 

Theorem 5.1. Suppose that two initial measures fiQ, uq € A^([0,27r) x M) and K satisfy 



(i) < De{i^o) < Deifio) <7T, / 9 i^io{de , dn) = 9uo{d9,dn) = tt. 

'[0,27r]xM J[0,2it]xR 



(ii) K > £'n(^to)max| . ^ — r, ■ }. , . \ 
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and let Ht and vt be two measure valued solutions to (jl.ip - (|1.2p corresponding to initial 
data fiQ and vq, respectively. Then, there exists to > such that 
r 2K cos D°° 1 

Wp(/X(,i/t) < exp {t-to) Wp{iJ.to,i^to), t>to, l<p<oo. 

L TT J 

Proof. First, we consider the case where p € [IjCxo). Note that the Wasserstein distance 
in one-space dimension is equivalent to the L^-distance of its corresponding pseudo inverse 
distribution function. Remember that we are assuming that the solutions are smooth, hence 
it is more convenient to obtain the L^-estimate from equation ()5.3p . Denoting by ^j, i = 1, 2 
the pseudo inverse functions associated to fit and i^t respectively, we get 

dt4>i = Q, + K / s'm{4>i^ — (j)i)dr]^dQ,^: , 







for i = 1,2. Then the above equations imply that 
dtih - (1)2) = K ( sin(</)i* - </)i) - sin((/)2* - (j)2))dr]^d9.^, 







IK I I COS y ^ 1 ^ ) sm ( ) drj^dih 



(5.6) 



^^\m-92UL. 



2pK 




, ,'i*-(pi , (P2*-(P2\ . 91* -91 92* 
COS 1 sm 



We multiply (15. Gh by psgn((/)i — (t)2)\<t'i ~ 4'2\^ ^ and integrate over [0,g{Q)] x M using the 
symmetry (r/, $7) <^=^ {r]^,Q^) to obtain 

d IP 

p 

rg{n) i-g{n,) 

COS I , , , 

2 2 J \ 2 2 

X [\(f>i - (f)2\P~^sgn{(f>i - (f)2) - \(f>i* - (f>2*\''~^sgn{(f)u - (t>2*)] drj^^drjdVl^^dVL. 
It follows from the proof of Lemma 15.41 that for all a, 6 G M, 

{\a\P-^sgn{a) - \h\P-^sg-a{b)) sin < 0. 

On the other hand. Lemma 1 5 . 1 1 implies that there exists to such that 

Dein) < D^, Deiut) <D^, t> to, 
and we use Lemma 15.21 to obtain 



max (j)i{g{n),n,t) - min (j)i(0,n,t) < D°° , 
nespt(g) n6Spt(g) 

max (j)2{g{n),n,t) - min (j)2{0,^,t) < D°° , t > to. 
r2espt(g) Qespt(3) 



Then, this yields 



n ^ noo ^ , r>l* - 91 , 92* 
< cos D < cos 1 



Hence, we obtain 

j^W'Pi - MIp < 2pKcosD°°J 

where 

J := / / sin 

JrxKjo jo 



>1* - 92* 
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X [I'/'i - <p2f ^sgn((/)i - 4)2) - - 4>2*\^ ^sgn - 4)2*)] dr]*dr]dil*dn. 
If we set <^ := — 02, then 

1^* - ^1 < |<^i* -4>i\ + 102* - 02| < 2D°° < vr, t > to- 
Since ^o, have the same center of mass, it fohows from Lemma 15.31 that 

^dridfl = 0, t > to- 







Thus, we can apply Lemma 15.41 with <I> = (/>i — 02 to obtain 
-^>„.,)< 



T.W^it^t, ^t) < W!^{^iu yt), t > to- 



This yields 

~ / 2K cos D°° \ ~ 

Wp{^it,i^t) <exp(^ {t-to)jWp{fito,i^to)- (5.7) 

In the case of p = cx), we use (j5.4p and (j5.7p to obtain 

~ / cos \ ~ 

Woo (/^t , i^t) <exjp\^ (t - to) j Woo ipto ,t^to)- 

This completes the proof for smooth solutions. As mentioned above a simple approximation 
argument as in Subsection 3.3 finishes the proof for measure valued solutions. □ 

Remark 5.1. The assumption in Theorem \5.1\ on the initial measures to have equal mean 
in 9 is not restricted. Due to Lemma \5.3l the mean in 6 is preserved in time. Thus, we can 
always restrict to the equal mean in 6 case by translational invariance of (jl.ip . 
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